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This  report  develops  a  detailed  analysis  of  a  type  of  electrohydro- 
dynamic power  generator  which  employs  an  ejector  and  a  so-called 
"fluid  flywheel"  as  essential  components.   The  medium  is  steam 
containing  electrically  charged  water  droplets . 

The  analysis  takes  into  account  the  experimentally  established  facts 
that  the  maximum  strength  of  the  electrical  field  that  can  be  sus- 
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tained  at  incipient  breakdown  at  the  most  critical  location  is  pro- 
portional to  the  fluid  density  at  that  location.   It  is  shown  that 
as  a  consequence  of  this  fact,  the  electrical  output  can  be  maxi- 
mized by  designing  the  primary  jet  for  an  exit  Mach  number  of  0.71. 

Estimates  are  made  of  the  pump  work  required,  of  mixing  losses  in 
the  ejector  and  of  friction  and  secondary  flow  losses.   The  mathe- 
matical analysis  is  reduced  to  a  fully  non-dimensional  form  and 
the  key  dimensionless  parameters  that  govern  performance  are 
clearly  identified.   A  preliminary  estimate  is  made  of  the  numerical 
values  of  these  parameters  and  the  overall  performance  of  the  system 
is  estimated  on  this  basis. 

Unfortunately,  the  results  so  obtained  are  very  pessimistic.  They 
indicate  that  even  at  100  atmospheres  pressure,  electrical  break- 
down so  severely  limits  power  output,  that  it  is  probably  insuffi- 
cient to  cover  the  demand  for  pump  power  and  for  power  to  overcome 
the  various  losses  involved,  let  alone  provide  any  useful  net  out- 
put. The  only  hopeful  note  is  that  the  analysis  so  clearly  pin- 
points the  problem  that  it  might  in  the  end  also  suggest  the  means 
for  surmounting  it. 

A  modified  design  where  the  electrical  conversion  section  is  placed 
after  the  ejector  is  analyzed  in  the  Appendix.   This  modification 
yields  less  pessimistic  results,  but  still  leaves  doubt  as  to  the 
viability  of  this  type  of  power  generator. 
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ABSTRACT 

This  report  develops  a  detailed  analysis  of  a  type  of 
electrohydrodynamic  power  generator  which  employs  an  ejector 
and  a  so-called  "fluid  flywheel"  as  essential  components. 
The  medium  is  steam  containing  electrically  charged  water 
droplets. 

The  analysis  takes  into  account  the  experimentally  estab- 
lished fact  that  the  maximum  strength  of  the  electrical  field 
that  can  be  sustained  at  incipient  breakdown  at  the  most  criti- 
cal location  is  proportional  to  the  fluid  density  at  that  loca- 
tion.  It  is  shown  that  as  a  consequence  of  this  fact,  the 
electrical  output  can  be  maximized  by  designing  the  primary 
jet  for  an  exit  Mach  number  of  0.71. 

Estimates  are  made  of  the  pump  work  required,  of  mixing 
losses  in  the  ejector  and  of  friction  and  secondary  flow 
losses.   The  mathematical  analysis  is  reduced  to  a  fully  non- 
dimensional  form  and  the  key  dimensionless  parameters  that 
govern  performance  are  clearly  identified.   A  preliminary 
estimate  is  made  of  the  numerical  values  of  these  parameters 
and  the  overall  performance  of  the  system  is  estimated  on  this 
basis . 

Unfortunately,  the  results  so  obtained  are  very  pessimistic 
They  indicate  that  even  at  100  atmospheres  pressure,  electrical 
breakdown  so  severely  limits  power  output,  that  it  is  probably 
insufficient  to  cover  the  demand  for  pump  power  and  for  power 
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to  overcome  the  various  losses  involved,  let  alone  provide  any 
useful  net  output.  The  only  hopeful  note  is  that  the  analysis 
so  clearly  pinpoints  the  problem  that  it  might  in  the  end  also 
suggest  the  means  for  surmounting  it. 

A  modified  design  where  the  electrical  conversion  section 
is  placed  after  the  ejector  is  analysed  in  the  Appendix.   This 
modification  yields  less  pessimistic  results,  but  still  leaves 
doubt  as  to  the  viability  of  this  type  of  power  generator. 
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1.   Introduction 

The  purpose  of  this  analytical  report  is  to  develop  a  mathe- 
matical model  of  a  proposed  type  of  electrohydrodynamic  (EHD)  power 
generator  which  embodies  an  ejector  as  an  essential  component.   The 
overall  configuration  is  approximately  toroidal  in  form  in  order 
to  help  conserve  the  angular  momentum  of  the  working  fluid  which 
circulates  through  it.   For  this  reason  the  inventor*  of  this  con- 
figuration has  dubbed  it  a  "fluid  flywheel". 

Fig.  1.1  is  a  schematic  diagram  of  the  device.   Superheated 
steam  at  stagnation  pressure   p   and  stagnation  temperature  T 
is  supplied  to  a  nozzle  which  discharges  into  the  ejector  mixing 
region  at  station  1.   A  secondary  fluid  enters  the  mixing  region  at 
station  2.   The  resulting  mixture  leaves  the  ejector  at  station  3. 
The  ejector  action  creates  a  pressure  rise  across  the  mixing  sec- 
tion which  is  just  sufficient  to  offset  the  effect  of  the  power  out- 
put and  of  the  various  flow  losses  that  occur  in  the  rest  of  the 
circuit. 

The  mixture  of  liquid  and  vapor  which  leaves  station  3  is  cooled 
during  its  return  circuit  sufficiently  to  condense  and  separate  out 
a  mass  flow  of  liquid  equal  to  the  mass  flow  of  steam  which  passes 

through  the  primary  jet. 

The  liquid  withdrawn  from  the  condenser  is  pumped  back  into  the 
boiler  and  superheater  where  it  is  heated  back  to  the  initial  stag- 
nation state  0 .   The  primary  fluid  thus  undergoes  a  Rankine  cycle 
where  the  turbine  and  condenser  have  been  replaced  by  an  EHD 
generator. 

*  See  note  in  Section  9,  References. 
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Conditions  in  the  primary  nozzle  are  controlled  in  such  a 
way  as  to  cause  fine  droplets  of  moisture  to  condense  here  and 
electrodes  are  provided  which  electrically  charge  the  liquid  drop- 
lets.  These  electric  charges  are  removed  from  the  flow  by  a  grid 
type  collector  located  at  the  ejector  exit,  station  3.   The  moving 
charges  constitute  an  electric  current  which  sets  up  a  correspond- 
ing electric  field.   The  electric  forces  acting  on  the  charged 
particles  tend  to  oppose  the  motion.   However,  the  size  and  charge 
of  the  particles  is  controlled  in  such  a  way  that  they  have  negli- 
gible mobility  with  respect  to  the  surrounding  fluid.   Thus  the 
moving  gas  does  work  on  the  fluid  particles  in  moving  them  down- 
stream against  the  resistance  of  the  electrical  forces.   In  per- 
forming this  work,  the  gas  stream  undergoes  a  corresponding  de- 
crease in  enthalpy.   The  work  done  by  the  gas  against  the  elec- 
trical forces  sets  up  an  electrical  potential  difference  between 
the  charging  electrode  and  the  collector  grid.   If  these  two 
terminals  be  connected  by  an  external  electrical  circuit,  that 
circuit  will  therefore  deliver  useful  electrical  power.   The 
only  turbo  device  required  is  a  small  pump. 

The  question  that  arises  is  whether  such  a  hypothetical  scheme 
can  be  made  practically  effective.   There  are  various  formidable 
obstacles  to  overcome,  including  the  achievement  of  negligible 
mobility.   It  is  therefore  important  to  have  a  complete  quan- 
titative model  of  all  cycle  components  to  help  guide  research 
and  development  efforts. 
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Such  a  model  has  now  been  developed  for  the  fluid  flywheel 
concept  and  is  presented  in  this  report.   The  model  shows  that  the 
performance  that  can  be  attained  is  governed  by  just  a  few  basic 
dimensionless  parameters.   It  identifies  these  parameters  and 
provides  some  basis  for  estimating  their  probable  magnitudes. 
By  identifying  the  key  parameters  and  showing  their  effect  on 
overall  performance,  this  analysis  permits  current  experimental 
research  efforts  to  be  redirected  along  more  fruitful  lines. 

This  report  is  the  result  of  research  sponsored  by  the  Depart- 
ment of  Energy. 
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2.   Basic  Electrical  Performance 

It  is  convenient  in  the  present  analysis  to  introduce  a 
parameter-  a      which  is  used  to  distinguish  between  the  two  dis- 
tinct cases  of  positively  and  negatively  charged  particles. 
Specifically,  we  set  a    =  +1   for  the  case  of  positively  charged 
particles,  and  we  set  a    =  -1   for  the  case  of  negatively  charged 
particles.   It  is  then  appropriate  to  denote  the  electrical  charge 
per  unit  mass  by  the  product  aq  where   q   has  the  units  of 
coulombs/kg  and  is  always  positive  by  definition. 

The  electrical  working  section  of  the  generator  coincides 
with  the  mixing  region  of  the  ejector  and,  assuming  polar  sym- 
metry, the  distribution  of  q   is  essentially  two  dimensional 
over  this  region.   However,  our  purpose  in  the  present  section  is 
to  develop  certain  key  dimensionless  parameters  which  depend 
primarily  on  how  q   varies  with  respect  to  the  axial  coordinate 
z.   These  parameters  are  not  sensitive  to  variation  with  respect 
to  the  radial  coordinate   r.   Under  these  circumstances  it  is 
permissible  to  simplify  the  analysis  in  the  following  way:   At 
each  axial  station   z   that  lies  within  the  mixing  region,  we 
subdivide  the  total  cross-sectional  area  A   into  an  inner  circular 
region  of  area  A'   and  an  outer  concentric  annulus  of  area  A" 
such  that 

A1  +  A"  =  A  =  A.  =  constant  along  channel       (2.1) 

While  fluid  density   p   (kg/m)  is  treated  as  essentially  uniform 
over  the  entire  mixing  region,  the  axial  velocity  (m/sec)  is 
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assumed  to  have  the  value   V   over  the  inner  area   A'   and  the 
value   V"  ■  Vp   over  the  outer  annular  area   A".   Similarly,  charge 
density  (coulomb/kg)  is  assumed  to  have  the  value   q'   over  the 
inner  area   A'   and  the  value   q"  =  0   over  the  outer  annular 
area   A".   Moreover,  the  quantities   A',  V  and  q1   are  regarded 
as  functions  of  coordinate   z   which  vary  monotonically  from  their 
initial  values   A, ,  V  ,  q    at   z  =  0   to  their  final  values 
A-,  V  ,  q    at  z  =  1 .   The  exact  form  of  these  functions  is  not 
of  immediate  importance  in  the  present  development. 

It  is  useful  to  introduce  the  following  auxiliary  definitions 
at  this  point  on 

I   =  length  of  working  section,  m 

V  =  magnitude  of  voltage  change  across  working  section, 
positive  by  definition,  volts 

z   =  £   =  permittivity  of  free  space 

=   8.854  x  10'12  farad/m 

<J>  =  electric  potential,  volts 

f  a  I  a  dimensionless  electric  potential 

S  =  y  =  dimensionless  axial  coordinate 

q1 
g  =  -r  =  dimensionless  charge/mass  ratio 
ql 

Poisson's  equation  and  its  first  two  integrals  can  now  be 
represented  in  the  following  one  dimensional  version: 
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c 


0     '     , c  c 

(♦-♦.)  =v(f-f0)  =  (||)^-!!i!iiy|g(c)d5d5  (2.4) 

0  0 


In  these  equations,  the  quantity   JL  is  the  negative  of  the 

dz 

electrical  field  strength  (volts/m) .   This  quantity  attains  its 
maximum  magnitude  at  station   z  =  0,   which  is  the  critical  loca- 
tion at  which  electrical  breakdown  first  occurs. 

It  has  been  established  by  experiment,  Ref .  (1) ,  that  over 
a  certain  range  of  pressures,  the  breakdown  strength  of  a  given 
medium  is  well  approximated  by  the  following  expression 


ldz/o 


Co  +  CBRpl 


(2.5) 


where   R   (joules/kg  k)  1S  the   gas  constant  of  the  medium  and  where 

2 

the   breakdown    constants      C       (volts/m)       and      C_.    (m  °K/coulomb)    are 

O  o 

characteristic  properties  of  the  medium  which  can  be  measured  ex- 
perimentally . 

The  evaluation  of  Eqs  (2.3)  and  (2.4)  at  the  exit  station 
z,    =  1   is  of  particular  importance.   In  this  connection  it  is  conven- 
ient to  introduce  the  following  notation: 


1 
iJj1  =  J   g(C)d5  (2-6) 
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i  e 

Hj2    =  /  J   gU)d*;dC 
0  0 


(2.7) 


Notice  that   <K   and   'K   are  dimensionless  coefficients  that 
depend  only  on  the  form  of  the  dimensionless  charge/mass  function 

g(C)  - 

The  length  I     of  the  working  section  of  an  ideal  EHD  generator 
having  a  given  current  flow   i   (amps)  should  be  so  chosen  as  to 
maximum  the  potential  difference   V   (volts)  and  hence  the  electri- 
cal power   P    (watts)   which  is  generated.   This  means  that  the 
e 

dimensionless  potential   f   should  attain  its  maximum  value  at 
the  exit  station,  or  that 


df 

dC 


=  0   at   5  =  1 


(2.8) 


It  also  follows  that 


(f  -  f J  =  a      at   c  =  1 


(2.9) 


Upon  substituting  Eqs  (2.5)  through  (2.9)  into  Eqs  (2.3)  and 
(2.4)  we  can  readily  extract  a  pair  of  useful  relations  from  the 
result.   In  this  connection  it  is  also  convenient  to  introduce 
the  auxiliary  notation 


a  = 


co/T. 


vCBR0l. 


\ 


2   Ml 


1 

TFTT 


(2.10) 
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Note  that  a   is  dimensionless .   This  quantity  is  normally  very  small 
compared  with  unity. 

Using  this  notation  and  procedure,  we  may  now  summarize  Eqs 
(2.3)  and  (2.4)  in  the  form 


(1+g)  (2.11) 

ql*      ty         £LBK 


/     ^2\  (2  12) 

v  =  V  "  *TJ(1+a)cBRpi" 


The  current  flow   i  (amps)   through  the  channel  may  be  written 

i  =  aq1P1V1A1  (2.13) 

The  gross  electrical  power  output   P   (watts)   is  now  given 
by  the  simple  expression 

P   =  [aV]i  (2.14) 

e 

In  order  to  reduce  our  final  result  into  a  more  useful  form 
we  also  write 


V1  =  a]M  (2.15) 

where 
V   =  axial  velocity  to  nozzle  exit,  m/sec 

a.  =  sonic  velocity  at  nozzle  exit,  m/sec 

NL  =  Mach  number  at  nozzle  exit 
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Upon  combining  the  last  five  relations  and  rearranging,  we 
obtain  the  important  result 

Pe  "  Ixf  "  ^)'1+a>2eCBR2eK«iAl  (2-16> 

For  the  purpose  of  non-dimensionalizing  this  result,  it  is 
useful  to  introduce  the  following  additional  dimensional  reference 
parameters 

a   =  sonic  velocity  corresponding  to  stagnation  conditions 
at  inlet  to  primary  jet,  m/sec 

P   =  density  corresponding  to  the  above  stagnation  conditions, 
kg/m3 

A   =  throat  area  of  primary  jet,  m 

2  2  2 

Now  dividing  Eq  (2.16)  through  by   eCgR  p0aQAt 

gives 


■  -  Wi  -  uUi^VteWV' 


e 


^R2povt/  hk  hi  r  '  \~o/  w^t 


Let  us  denote  the  bracketed  quantity  on  the  right  in  the 
following  way 

;(1+a)2(^)2(^)Mi(^t) r  a+a)2]p3  f(mi}  (2-18) 

where   ij>3   is  a  certain  normalizing  constant  which  will  be  defined 
presently. 
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In   order   to    analyze   the    function   defined   by   Eq    (2.18),    we   treat 
the   expansion   through    the   primary   nozzle   in   the    usual  way   as    an 
isentropic   process    of   constant   total    enthalpy.      The  medium  is 
treated   as    a   perfect  gas    for  which      y    ,    the   ratio   of   specific  heats, 
is    a   constant. 

The    relations   which   govern   such   a   nozzle    flow   are  well   known 
and  will    therefore  be   used    freely  here  without   detailed  deriva- 
tion.     Note   that   there    are    two   distinct   expressions    for    the    func- 
tion     i|;3F(M,)       depending  on  whether     M,       is    subsonic   or   supersonic. 
For   the   subsonic   case,    the    throat  of    the   nozzle   coincides   with    the 
exit   station   so    that      (A, /A.)       equals    unity.      In   the   supersonic 
case       (A- /A.)      becomes    a    function   of      M, . 

When   these   details    are   properly  worked  out,    the    following 
solution   is    obtained. 


For 


M,   <    1 


i|»3F(M1)     =   Mx 


1*151 


M- 


(V  +  3) 

2(y-i: 


(2.19 


For 


M,    >    1 


i^F(M. 


(y+1) 

JTFI1    r 


Y  +  1 


1    + 


^K 


(Y-D 


(2.20) 
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By  differentiating  Eq  (2.19),  it  can  be  shown  that  the  function 
ijj^F  (M^  )   passes  through  its  maximum  value  precisely  at 


Mlopt  =  J  *  °-71  (2.21) 


Note,  incidentally,  that  this  result  happens  to  be  independent 

of  y[ 

2 

On  the  other  hand  the  function   (1  +  a)   <Jj_F(M,  )   does  not 

necessarily  pass  through  its  maximum  at  precisely  this  same  Mach 

number  because,  as  may  be  seen  from  Eq  (2.10),  the  quantity   a 

is  itself  a  function  of  density   p   at  station  1,  and  hence  of 

M..  .   Nevertheless,   a   is  very  small  compared  with  unity,  and  its 

variation  will  therefore  have  an  effect  on  Mn      so  small  ^hat 

lopt 

it  can  safely  be  neglected. 

As  a  normalizing  condition  we  stipulate  that  the  maximum 
value  of  the  function  F(M.)  shall  be  unity.   Thus 


F     =  F(M.   J  =  1  (2.22) 

max       lopt  \*.**i 


This  constraint  now  fixes  the  magnitude  of  the  normalizing  con- 
stant  if>~  which  works  out  to  be 


(Y+3) 

"           -l' 

1 

4 

*3 

V2~ 

Y+3 

(2.23) 
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For  example,  for   y  =  1-3   this  formula  gives 

^3  =  0.4211  (2.24) 

The  dimensionless  coefficients   tK,  ^2  and  ty  2      and  the  factor 

2 
(1  +  a)   can  next  be  combined  into  an  overall  dimensionless  coef- 
ficient  tp  ^   according  to  the  relation 

^3  /    *2\ 
U    =~  (l-^jU+a)2  (2.25) 

With  this  notation,  the  electrical  power  output  can  be  expres- 
sed in  the  following  form,  namely, 

\ 

2  2   2 1"  *4F(M1» 

where  the  function   F(M  )   is  defined  by  Eqs  (2.19),  (2.20) 
and  (2.23) . 

The  function   F(M  )   is  listed  in  Table  2.1  and  plotted  in 
Fig  2.1.   These  results  are  for   y  =  1.3  which  is  a  typical  value 
often  used  for  steam.   Notice  that   F(IXL)  reaches  its  peak  value 
of  unity  at  an  optimum  exit  Mach  number  of  approximately   0.71. 
Observe  that  appreciable  deviations  from  this  optimum  exit  condi- 
tion cause  heavy  performance  losses . 

Under  these  conditions,  the  subsequent  analysis  will  be  re- 
stricted specifically  to  the  optimum  exit  Mach  number  M     =0.71 
as  derived  above.   Consequently,  Eq  (2.25)  can  then  be  simplified 
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TABLE  2 . 1 
y  =  1.3 


M-l     0.0  F(MX)     0.0 

0.20  ====         0.4550 

0.40  0.8014 

0.60  0.9774 

0.71  1.0000 

0.80  0.9849 

1.00  0.8722 

1.20  0.7241 

1.40  0.5886 

1.60  0.4704 

1.80  0.3711 

2.00  0.2901 

2.20  0.2253 

2.40  0.1744 

2.60  0.1347 

2.80  0.1041 

3.00  0.0805 

3.20  0.0625 

3.40  0.0486 

3.60  0.0380 

3.80  0.0298 

4.00  0.0235 

4.20  0.0186 

4.40  0.0148 

4.60  0.0119 

4.80  0.0095 

5.00  0.0077 
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Fig.    2.1      F(M1)       versus      M 
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to  the  form 


Z  lS    2 ~|~  ]h  (2.27) 

eCn  R  Pn  a^A. 


In  order  to  get  the  greatest  benefit  from  Eq  (2.27)  we  should 
know  the  value  of  the  power  coefficient  ty.   .   This  coefficient 
is  defined  by  Eqs  (2.6),  (2.7),  (2.23)  and  (2.25).   It  may  be  seen 
from  Eqs  (2.6)  and  (2.7)  that  in  order  to  calculate  the  quantities 
'^,   and  i)~      that  ultimately  fix  the  initially  unknown  coefficient 

\\).,   we  must  know  the  exact  dimensionless  charge/mass  distribution 
as  fixed  by  the  function   g(?) .   Unfortunately  this  function  re- 
mains unknown  so  that  an  exact  solution  is  not  possible  within  the 
scope  of  the  present  simplified  analysis. 

Nevertheless ,  we  do  have  some  approximate  information  about 
this  function  which  is  adequate  to  provide  a  first  order  estimate 
of   'J;..   Let  us  call  this   ^4n  •   In  order  to  establish  such  an 

estimate  or  reference  value  it  is  only  necessary  to  assume  a 
qualitatively  correct  and  analytically  convenient  form  for  this 
function. 

Note  that 

g  =  1  at  c  =  0  (2.28) 

_  ^3 

9  '  g3    ql  at  r,    =  1  (2.29) 

Moreover   g(C)   varies  monotonically  between  these  two  end 
values  and  it  should  also  satisfy  the  condition,  at  least 
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approximately,  that 


if  .  0   .t  c  -  1 


It  is  also  timely  to  note  here  that  the  current  i   evaluated 
at  the  two  stations  1  and  3  gives  the  relation 

qlplAlVl  =  q3p3A3V3  (2.31) 

From  this  we  readily  infer  that 


q3        /PjAjVA 

93  =  ~i  mvm " x  (2-32) 


where   x   is  the  dimensionless  mass  flow  ratio  of  the  ejector 
A  simple  analytical  form  for  the  reference  function   g(C) 
that  satisfies  the  foregoing  constraints  is 


g(C)  =  x  +  (i-x)(l-c) 


2  (2.33) 


Upon  substituting  this  into  Eqs  (2.6)  and  (2.7)  and  inte- 
grating we  obtain  the  expressions 

'h    =   1    (1+2x)  (2.34) 


ty2    =  j    (1+x)  (2.35) 
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Then  using  Eq  (2.25)  we  obtain  the  desired  reference  value  as 


40 


j   *3(l+o) 


(l+5x) 
(l+2x) 


(2.36) 


where  constants   ^_   and   a   are  still  defined  by  Eqs  (2.23)  and 
(2.10)  as  before. 

Of  course  ty.Q      is   only  a  reference,  not  the  true  value  \\>. 
itself.   However,  we  may  relate  these  two  quantities  as  follows. 


]h    =    u^40 


=  U 


4  *3 


1+ct) 


(l+5x) 
J  (l+2x) 


(2.37) 


Here   y   is  regarded  as  an  empirical  correction  factor  of  order 
unity.   It  is  to  be  expected  that  whereas   i|»4  itself  may  tend  to 
vary  over  a  certain  range,   y   should  be  more  nearly  constant. 
We  may  now  combine  Eqs  (2.27)  and  (2.37)  in  the  form 


2   2 
yec  *ITp   a  A., 
B     o   o  t 


|  <Ml+a)2 


(l+5x) 
(l+2x) 


(2.38) 


The  quantities   ^3,  a  and  x  which  appear  in  Eq  (2.38)  are 
defined  by  Eqs  (2.23),  (2.10)  and  (2.32),  respectively.   The  coef- 
ficient  u   must  be  determined  by  experiment  or  by  more  advanced 
analysis.   However,  for  purposes  of  order  of  magnitude  estimates 
it  is  permissible  to  take 


y  -  1 


(2.39) 
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3.   Gross  Electrical  Power:   Change  of  Dimensional  Base 

In  analyzing  the  mixing  and  friction  losses  in  the  ejector 
and  the  friction  and  secondary  losses  in  other  parts  of  the  sys- 
tem, it  will  prove  useful  to  employ  as  dimensional  reference 
quantities  the  fluid  density   p=p-,=p2  =  P3   in  tne  ejector, 
the  exit  area   A.,   from  the  ejector  and  the  average  axial 
velocity  V,   at  this  station.   Thus  the  reference  mass  flow  rate 

becomes   p.A.V   (kg/sec),  the  reference  kinetic  energy  per  unit 

2 
mass  becomes   V_/2   (joule/kg)  and  the  reference  power  becomes 

2 
(p-A-Vj  (V,/2)  (watts).   In  this  connection  it  is  appropriate  to 

define  a  dimensionless  gross  electrical  power  coefficient   c 

in  the  following  way,  namely, 


P 

c   =  — _!   .  -  (3.1) 

e 


[P3A3V3][V32/21 


It  is  also  necessary  to  establish  suitable  links  between  the 
above  reference  quantities  and  certain  corresponding  quantities 
at  station  1,  the  exit  from  the  primary  jet.   In  this  connection 
recall  that 


m,\  /piA,V,\ 

— ±.\=   •*■  L    x  =  v  =  dimensionless  mass  flow  ratio        (3.2) 

v*3/  \p3A3V3/ 
We  also  elect  to  introduce  the  following  definition 
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'Vl\  (3.3) 

tt~ )-   Y  =  dimensionless  velocity  ratio 

VV3/ 

As  before  it  will  be  convenient  to  write 


Vl  =  alMl  (3*4) 


Moreover  we  can  substitute  for  P  in  Eq  (3.1)  its  expression 
as  previously  developed  in  Eq  (2.38). 

When  the  foregoing  relations  are  combined  and  rearranged,  the 
following  result  is  obtained 

|     —    (1+cO2  MyeC^P^p    /  k,,^       , 

)  4      Y (JT       B         Ho(Ml+5x)xy    y  ,,    cn 


This  suggest  the  utility  of  defining  the  following  auxiliary 
quantities,  namely, 

if^«>2       ) 


'KSfe)Ml3fe)l  '3-6» 


and 

2  2 

^£C3  R  po  x 

2 1 

a 
o 

Recall  that  the  quantities  involved  in  Eq  (3.6)  all  correspond 
to  the  value   M]_opt  =  l/v'2.   Hence  ^       turns  out  to  be  a  constant 
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which  depends  only  on  y   and  a.   When  the  details  are  worked  out 
we  find  that 

/  y  +  i\ 


*  =  ir   4 


Y  Y  +  3 


Y  "  1 


d  +  <*)2  (3.8) 


For  example,  for  y   =1.3   this  gives 

^  =  1.32(1  +  a) 2  (3>9) 

Since  the  medium  is  treated  as  a  perfect  gas,  the  crucially 
important  parameter  defined  by  Eq  (3.7)  can  be  further  reduced  to 
the  form 

ZBjc 

T  2 
0  (3.10) 

This  result  shows  the  importance  of  maintaining  inlet  stagna- 
tion pressure   p    as  high  as  possible  and  inlet  stagnation  tempera- 
ture  T    as  low  as  possible.   However,  if  condensation  is  used  to 
o  c 

produce  droplets,   T    is  restricted  to  that  range  of  temperatures 
which  produces  optimum  droplet  size.   One  way  around  this  restriction 
is  as  follows:   Let  the  carrier  fluid  be  changed  from  steam  to  air 
(or  to  some  other  gas  which  is  noncondensible  at  ordinary  tempera- 
tures) .   Let  charged  liquid  droplets  or  solid  particles  of  proper 

size  be  injected  into  the  carrier  fluid.   This  permits   T    to  be 
J  c  o 

sharply  reduced  whereupon   S   is  correspondingly  increased.   Of 
course  this  alternative  may  involve  its  own  characteristic  draw- 
backs as  well,  but  these  will  not  be  discussed  further  at  this 
point. 
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While  the  introduction  of  the  separate  quantities   y,  if;  and  6 
follows  naturally  from  Eq  (3.5),  we  note  that  these  quantities 
always  occur  in  the  form  of  the  product  y^g  .   Hence  it  will 
simplify  the  subsequent  development  to  set 


u  =  mp3 
whereupon  Eq  (3.5)  reduces  to 


(3.11) 


n   <1+5X)       2 

ce  =  u  2  *y 


(3.12) 


(1+2X) 

In  the  subsequent  analysis  of  ejector  performance,  it  will  be 
shown  how  the  two  fundamental  parameters   x  and  y   govern  pressure 
rise  through  the  ejector,  mixing  losses  and  so  on.   Eq  (3.12) 
therefore  ties  in  gross  electrical  power  output  with  these  other 
effects  on  the  same  consistent  overall  basis.   In  particular,  it 
is  important  that  the  gross  electrical  power  coefficient  ce    as 

expressed  by  this  equation  should  exceed  the  sum  of  all  the 
losses  in  the  system  by  as  wide  a  margin  as  possible. 
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4,   Ejector  Performance 

Recall  the  notation  introduced  earlier,  namely, 


-.  (4.1) 

mass  flow  ratio 


lv~~/  =  y  =  velocity  ratio 

As  a  basis  of  comparison  with  the  actual  physical  ejector, 
consider  a  hypothetical  ideal  ejector  which  satisfies  the  following 
constraints.   Firstly,  the  fluid  density  is  uniform  throughout 
the  ejector  so  that   p-,  =  p2  =  C3  =  o  =  constant.   Secondly,  pres- 
sure is  uniform  across  the  entire  inlet  so  that  p_  =  p . .   Thirdly, 
the  velocities  are  purely  axial  and  are  uniform  across  sections 
1,  2  and  3.   Fourthly,  wall  friction  is  negligible. 

Imagine  this  ideal  ejector  to  be  in  operation  with  all  elec- 
trical circuits  turned  off.   The  ideal  pressure   (p-,*  -  p-,)   that 
would  occur  under  these  conditions  may  be  expressed  by  means  of 
the  dimensionless  coefficient 


(p3*  "  P,)      *  (4.3) 

—= =  c    =  ideal  pressure  rise  coefficient 

ypV.,2        P     (for  given  values  of   x  and  y) 

Of  course  the  real  physical  unit  operating  under  realistic 
conditions  with  electrical  power  output  connected  to  a  load,  will 
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produce  some  smaller  pressure  rise   (p-,  -  p,  )   at  the  same  values 
of   x  and  y.   This  actual  pressure  rise  may  be  expressed  by  the 
coefficient 


3    ^1   =  c   =  actual  pressure  rise  coefficient        ,.    .* 
1  T7  2  (for  given  values  of  x   and  y) 

2pV3 


The  difference   (CD*  -  c  )  is  accounted  for  in  part  by  the  gross 
electrical  power  output  as  expressed  by  the  previously  analyzed 
coefficient   co   and  in  part  by  additional  mixing  and  friction 

losses  that  characterize  the  real  system. 

The  mass  flow  rates  across  sections  1,  2  and  3  may  be  written 

™1  =  pAlVl  =  xpA3V3  (4-5>  ■ 

m2  =  pA2V2  =  (1-x)  A3V3  (4.6) 

m3  =  pA3V3  (4-?) 

From  these  relations  and  from  Eqs  (4.1)  and  (4.2)  we  find 
the  two  area  ratios 


*3/     *  (4.8) 


rJ =  (1-x)fe)  (4-9) 
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These  ratios  satisfy  the  constant  area  condition 


a +  g  =  ^ 


We  can  solve  these  last  three  equations  for  the  velocity 
ratio 


V2\     (l-x)y 

vj     (y-x)  (4.1i: 


The  momentum  equation  can  now  be  applied  to  the  ideal  ejector 
as  follows 

(P3*  -  PX)A3  =  PA3V32   xDj  +  (l  -  X)(!ij  -  1  (4.12) 


2 
We  next  divide  through  by  ijPA-V-   and  use  Eqs  (4.2)  and  (4.11) 

to  eliminate  the  velocity  ratios  (V../V-).   and   ^V2//V3^  *   After 
simplification,  this  reduces  to  the  important  result 


c  *  =  lx(y  -  l)  (4.13 

P       (y-x) 
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Next  consider  a  purely  hypothetical  reversible  device  which 
is  not  an  ejector  at  all  but  a  type  of  ideal  turbomachine .   This 
device  receives  two  streams  of  fluid  identical  in  all  respects  to 
the  streams  at  stations  1  and  2  of  the  ideal  ejector.   It  dis- 
charges a  stream  identical  in  all  respects  to  the  stream  at  sta- 
tion 3  of  the  ideal  ejector  except  for  the  pressure   pR3   which 

* 
is  higher  than   p,   .   Again  the  density   p   is  taken  as  constant. 

The  reversible  device  operates  adiabatically  and  hence  isentropi- 

cally.   Under  these  conditions  the  following  energy  equation  holds 

Recall  that   p   =  p   . 

xm3(p1   +   I>dVi2)    +    (1   "   x)    ™3    (Pi   +   IpV22) 

=   S     (?3R   +   IpV32)  (4'14) 

The  aPove  equation  can  be  simplified  by  subtracting  from  it, 
term  by  term,  the  following  identity. 

xm3p1  +  (1  -  x)m3p1  =  ^3P1  (4.15) 

Upon  carrying  out  the  above  subtraction,  then  dividing  through 

2 
by  the  reference  power  m,  (V^  /2)   and  simplifying,  we  readily 

find  the  pressure  rise  coefficient  for  the  reversible  device  to 

be 


(P3R  "  ?3}      „  2    (1  -  x)V 


=   CpR  =    xy"  +  is -2   -  1  (4.16) 


IpV  2  /         {*   -   x) 

2P  3 
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Comparison  of  this  result  with  Eq.  (4.13)  shows  that  even  the 

ideal  ejector  has  a  certain  inherent  mixing  loss,  call  it   c 

m 

as  compared  with  a  truly  reversible  device.   Thus 
* 
°m   =  (CPR  ~  Cp*)  (4.17) 

Of  course  the  mixing  loss   c   of  the  actual  ejector  will 

m  J 

* 
be  somewhat  larger  than   c       , 

m   and  it  is  advantageous  to  express 

it  in  the  form 

m 
cm  "  n^~  (4.18) 

where   nE  ,    which  we  term  ejector  effectiveness,  is  a  dimension- 
less  parameter  smaller  than  unity  which  must  be  determined  from 
experimental  data. 

Eq.  (4.18)  identifies  one  of  the  component  losses  that  reduce 

the  actual  pressure  rise  coefficient   c   below  its  ideal  rever- 

*  P 

sible  value   c-,-  .   Another  loss  is  that  associated  with  Separa- 
te K 

tion,  friction  and  secondary  flow  effects  within  the  ejector;  let 
us  denote  it  by  symbol   c^E  .   We  have  seen  that  a  third  component 
is  simply  the  gross  electrical  power  output  as  expressed  by  the  co- 
efficient  c   .   Thus  the  overall  energy  balance  within  the  ejec- 
tor requires  that 

Cp  "  CPR  "  i  (CPR  "  Cp*}  '  CfE  "  Ce  (4'19) 

E 

The  resulting  pressure  rise  through  the  ejector  as  expressed 
by  Eq.  (4.19)  must  be  just  great  enough  to  offset  the  losses  in 
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the  return  passage.   These  again  are  losses  associated  with  separa- 
tion, friction  and  secondary  flows;  let  us  denote  them  by  coeffi- 
cient  cfR  .   Accordingly,  for  the  return  passage  we  may  write 


p     f  R 


(4.20) 


We  now  eliminate  c  between  Eqs .  (4.19)  and  (4.20)  and  bring 
all  terms  to  one  side.  We  also  introduce  the  following  definition 
for  the  total  friction  loss  in  the  fluid  flywheel,  namely, 


cf  ~  cfE  +  cfR 


(4.21) 


The  result  is 


inE  (CPR   CP 


+  c. 


-  c   =  0 

e 


(4.22) 


Now  using  Eqs.  (4.13),  (4.16)  and  (3.12)  to  express  the  above 
components  in  terms  of  variables   x   and  y  ,  we  obtain 


2 

xy   + 


m     i  3  2 

t1  -  x)  Y  -  i 

2 

(y  -  x)z 


E 


2 
xy  + 


(i  +  2Xr 


=  0 


(1  -  x)3y2    .    2x(y  -  l)2 

7~~T--  1 (y  -  x) 

\  y  -  x  / 


(4.23) 


This  is  the  overall  energy  balance  equation  for  the  "fluid 
flywheel"  and  it  is  fundamental.   Notice  that  the  term  in  the  first 
pair  of  curly  brackets  represents  the  ideal  power  output  of  a  re- 
versible device,  the  term  in  the  second  pair  of  curly  brackets 
represents  the  sum  of  the  fluid  power  losses  caused  by  mixing  and 
friction  and  the  term  in  the  third  pair  of  curly  brackets  represents 
the  gross  electrical  power  output. 
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Eq.  (4.23)  suggests  the  utility  of  defining  a  "conversion 
efficiency"   ricv  as  the  ratio  of  the  gross  electrical  power 
to  the  ideal  power   cpR  .   Thus 


c 


(1  +  5x) 


-2  xy 


ce  (1  +  2x) 

Ky2   +  (1  -  x)  J      -1 
(y  -  xr 

Notice  that  low  values  of  the  parameter   v  ,  which  is  ultimately 
determined  by  the  electrical  breakdown  limit,  will  necessarily 
produce  low  conversion  efficiencies. 

In  connection  with  Eq.  (4.23)  and  (4.24),  notice  that  x 
always  falls  between  zero  and  unity  while   y   is  always  greater 
than  unity  but  has  no  preassigned  upper  limit. 

Assuming  that  parameters   nE  /  cf   and  v   are  known  or  can 
be  estimated,  Eq.  (4.23)  then  fixes  a  corresponding  unique  relation 
between  variables  x   and  y  .   This  may  be  described  by  a  line  in 
the  xy   plane.   All  subsequent  calculations  must  be  confined  to 
points  that  lie  along  this  line. 

If  the  area  ratio   _1   be  already  fixed,  then  we  have  from 

Eq.  (4.8)  the  additional  constraint 
/A. 


A) =  & 


(4.25) 


Of  course  for  a  fixed  value  of   (A^/A^ ,  Eq.  (4.25)  also 
represents  a  line  in  the  xy   plane.   In  fact  it  is  a  straight 
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line  of  slope   (A3/A. )   which  passes  through  the  origin.   The  one 
or  more  points  of  intersection  of  this  straight  line  with  the 
curve  associated  with  Eq  (4.23)  fix  the  conditions  at  which  steady 
operation  of  the  fluid  flywheel  is  possible.   The  location  of  such 
a  point  determines  the  values  of   x   and  y   and  hence  ultimately 
determines  all  other  performance  characteristics  of  the  system  as 
well. 

Notice  that  Eq.  (4.2  3)  lends  itself  to  further  algebraic  sim- 
plification which  makes  the  subsequent  numerical  solution  easier. 
On  the  other  hand  the  terms  in  the  simplified  result  cannot  be 
given  the  clear  physical  interpretation  that  is  possible  with 
Eq.  (4.23)  itself.   After  simplification,  Eq.  (4.23)  becomes 

nE     (y-x)    ^   1)  |*    (y_  x)2 

(1  +  5x)     2         -.  i  a    -^c\ 

-   v xy   -  c-  =  0  (4.26) 

(1  +  2x) 

For  purposes  of  numerical  solution,  it  is  useful  to  cast 
Eq.  (4.26)  into  a  more  tractable  form  as  follows.   Let 

A  =  2(y-l)2 
l     nE  (4.27) 

*2m(k     ~    X)y2  (4'28) 


A.  =  vy2  (4.29) 
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A,  =  c~  -   =•     +   1 

4    f   nE 


(4.30) 


With  this  notation,  Eq.  (4.26)  may  be  rewritten  as  follows, 
where   F   is  a  function  that  must  equal  zero.   Thus 

3  i 


F  =  A 


X 


1  (y-x)   A2 


l-xi 


(y-x) 


-  A   Ii±5xjx  .  A   =  o 

3         2     4 

J  (l+2x)2     4 


(4.31) 


2         2 

Next  we  multiply  through  by   (1  +  2x)   (y-x)    and  denote 

the  resulting  function  by   G.   Upon  expanding  the  right  side, 
regrouping  terms  and  simplifying,  we  find  that   G   is  a  quartic 
in  x.   The  coefficients  of  this  polynomial  turn  out  to  be 


B„  = 


4  =  -  4A1-4(3-2y)A2-5A3-4A4 
33  =  4(y-l)A1-4y(y-2)A2+(10y-l)A3+4(2y-l)A4 


B„  = 


B„  = 


(4y-l)A1-(4y2-2y-5)A2-y(5y-2)A3-(4y2-8y+l)A4 


yA1-(y2+l)A2-y2A3-2y(2y-l)Ai 


Bn  =  -A?-y  A^ 


(4.32) 
(4.33) 
(4.34) 
(4.35) 
(4.36) 


The  function   G  may  then  be  written 


G  =  B4x4+  B3x3+B2x2+B1x^B0  =  0 


(4.37) 


The  unknown  root   x   can  be  found  by  Newton's  method.   For 
this  purpose  we  require  also  the  derivative  of   G  which  is 


G1  =  4B4x3+3B3x2+2B2x+B1 


(4.38) 
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Let  x    be  any  trial  value  of   x  and  let   G    and   G' 

n       2  n         n 

be  the  corresponding  values  computed  from  Eqs .  (4.37)  and  (4.38) 
Then  according  to  Newton's  method,  the  next  trial  value  x  +, 
should  be  taken  as 

xn+l=*n-^  (4-39) 

n 

This  procedure  converges  rapidly  to  the  true  root  x  . 
This  method  was  used  in  calculating  the  numerical  example  which 
is  presented  in  a  later  section. 
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5.   Thermodynamic  Performance 

As  shown  in  the  schematic  temperature  entropy  diagram,  Fig 
5;1,  the  tluid  flywheel  receives  a  supply  m   (kg/sec)  of  super- 
heated steam  at  condition   p   ,  T   and  discharges  saturated  or 
nearly  saturated  liquid  at  condition  p.  ,  T.  .   It  gives  up  heat 
to  the  ambient  atmosphere  at  absolute  temperature   T,   and  delivers 
a  gross  electrical  power  output  P    (watts) .   It  also  requires 
an  input  of  electric  power,  call  it  P   ,  to  energize  the  injec- 
tor.   This  is  normally  a  small  percentage  of  the  gross  electrical 
output  and  is  conveniently  expressed  by  means  of  an  excitation 

efficiency   n    such  that 
x 


P   -p=p»   =nP  (watts)  (5.1) 

e    x     e    'x  e 


^continued  on  page  5-3! 
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Fig.  5.2.   Schematic  Temperature  Entropy  Diagram  for  Primary 
Stream. 
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After  being  condensed  and  leaving  the  fluid  flywheel,  the 
liquid  water  must  be  pumped  back  into  the  boiler.   For  an  ideal 
pump  of  100%  efficiency  this  would  require  the  expenditure  of  shaft 
power  in  the  amount 

Pp   =  ™1  (po  "  p4]  v4         (watts)  (5.2) 

where   v4  denotes  the  specific  volume  of  the  liquid  at  state  4. 
For  a  real  pump  of  efficiency   n_   the  corresponding  power  is 

ir 
* 

P 

p  =  -£- 

P    np  (watts)  (5.3) 

Thus  it  is  seen  that  for  the  overall  cycle,  the  net  useful 
power  output  may  be  written  as 

V  =  <nxPe  -  pp*/np)        (watts)  (5.4) 

The  rate  of  heat  input  to  the  steam  passing  through  the  boiler 
and  superheater  is 

* 
Q  =  m,  (n  -h.)  -  P^  /n       (watts)  (5.5) 

1    O   4      P    P 

where   h    and  h.   denote  the  enthalpies  (joule/kg)  at  states 

O  4 

o  and  4 . 

For  the  hypothetical  case  of  100%  pump  efficiency,  we  denote  the 

•  * 
corresponding  heat  rate  calculated  from  Eq.  (5.5)  by  symbol   Q  . 

This  differs  only  very  slightly  from  Q  . 
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The  overall  cycle  efficiency  now  becomes 

n   =  -£-  =  i.(n  P   -  P  */n  J  (5.6) 

C    Q     Q  \  *  S     P    P/ 

It  is  of  interest  to  compare  the  above  performance  of  the 
actual  EGD  device  with  the  ideal  performance  of  a  hypothetical 
reversible  device  operating  under  the  same  input,  output  and  am- 
bient conditions.   According  to  the  theory  of  availability  in 
steady  flow,  the  maximum  electrical  power  that  could  be  produced 
by  such  an  ideal  device  is 

Pe*  =  ™1  )(ho  "  V  '  TA  (so  '  s4)l    (watts)  (5.7) 

where   s    and   s4   denote  the  entropies  at  states   o   and   4, 
respectively  (joule/kg  °K) . 

The  cycle  efficiency  of  this  ideal  device  would  be 

*    1     *  -  P  ) 
nc  -  t*  (P     P  ]  (5.8) 

Q 

It  is  also  instructive  to  compare  the  net  output  of  the  real 
system  with  that  of  the  above  ideal  and  reversible  device.   We 
term  this  ratio  the  relative  output.   Thus 

/n  P   -  P   /n 
n    -/  x  g    P/  P 

rel   \Pe*  "  Pp*   /  ^5-9) 

All  of  the  foregoing  relations  can  be  put  into  dimensionless 
form  by  dividing  all  quantities  having  the  units  of  power  by  a 
suitable  reference  power.   In  the  previous  analysis  of  the  fluid 
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flywheel,  it  was  natural  and  convenient  to  choose  this  reference 
as  the  kinetic  power  of  the  total  stream  at  the  ejector  exit, 
station  3.   While  it  would  be  possible  to  retain  this  same  reference 
in  this  section,  it  is  advantageous  to  choose  here  as  reference, 
the  kinetic  power  of  the  primary  jet  at  the  ejector  inlet,  station 
1.   This  choice  reduces  the  final  results  to  a  simpler  form  and 
makes  them  somewhat  easier  to  interpret.   Thus  the  reference  power 
is  now 

"A^y  »i  <ho  -  v  <5-10' 

In  computing  this  reference  power  it  is  permissible  to  assume 
that  state  1  is  fixed  by  the  two  conditions 


s,  -  s  (5.11) 

1    o 


P]_  =  p4  (5.12) 


In  connection  with  Eq.  (5.10)  it  is  also  instructive  to  note 
that 


V,2    \   a^2M,2 
(h   -  h,>  -  JL' 


1+I^V 


(5.13) 


Recall  also  that  if  we  introduce  the  optimizing  condition 
developed  earlier,  namely, 

M  2  =  H ,2         -  i  (5.14) 

1      lopt    2 
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then  Eq  (5.13)  simplifies  further  to 


V. 


(h  -  h.)  =  -t 

O     1      2 


(Y+3) 


(5.15) 


Incidentally,  it  is  also  worth  noting  that  the  overall  pres- 
sure ratio  across  the  fluid  flywheel  is 


11)1  11 


i+teiKa 


Y-l 


(5.16) 


Moreover,  if  we  observe  the  restriction  of  Eq  (5.16)  this  becomes 


P~ 


rp^> 

vp7, 


,Mr 


y  +  3 


(5.17) 


For   example,    for      y   -  1.3      this    gives 


vPj     \Po/~    °'731 


(5.18) 


This  result  expresses  one  of  the  fundamental  limitations  of 
the  fluid  flywheel  concept,  namely  the  restriction  to  a  fixed 
pressure  ratio.   The  reason  for  this  restriction  was  explained 
earlier  in  connection  with  Fig.  2.1. 
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Upon  simplifying  and  reducing  the  various  foregoing  thermo- 
dynamic relations  in  the  manner  explained  above,  we  obtain  the 

following  results. 
* 


m.  (h  -h, ) 
l  o   I 


*    )     (T4"TA) (so~s4) '   Reversible 

Ke    i        (h  -h.. )     I   Electric  Power 

o   1 


(5.19) 


m. (h  -h. ) 
1  o   1 


nu  (h  -h, ) 
1   o   1 


=  K 


(PQ-P4)v4 


=  K, 


<Vhi> 

(Vh4» 

* 

KP 

(h   -h.) 

Reversible 
Pump  Power 

Reversible 
Heat  Input 


5.20) 


(5.21) 


m1  (h  -h, ) 
1  o   l 


=  K, 


(hQ-h4>    ^ 


Actual  Heat 
Input 


(5.22) 


m  (h  -h. | 
1   o   1 


=  K   =   v 
e 


(l+5x) 
(l+2x) 


Gross  Electric 
Power 


(5.23) 


m.  (h  -h1 ) 

1   o   1 


=  K 


p 


Net  Electric 
Power 


(5.24) 


"°   =  4  VA    "  ^P~ 


Cycle  Efficiency   (5  2s) 


nc  =  A  (Ke  -  Kp 


Reversible 
Cycle  Efficiency 


(5.26) 


'rel 


,n  k  -  k  /n 
-  x  e    p   p 


Relative  Output    ,  5  .^-j) 
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Once  Eq.  (4.37)  of  the  previous  section  is  solved,  it  becomes 
a  relatively  straight  forward  matter  to  apply  the  relations  of 
the  present  section  to  calculate  the  various  thermodynamic  per- 
formance parameters.   Recall  that  Eq.  (4.37)  defines  a  whole  series 
of  points  in  the   xy   plane.   By  calculating  the  thermodynamic 
performance  parameters  of  this  section  for  such  a  series  of  points, 
and  by  comparing  the  respective  results,  it  becomes  possible  to 
choose  that  optimum  design  point  which  yields  the  best  overall 
trade  off  among  the  various  desired  performance  characteristics. 

So  far  only  an  initial  trial  calculation  of  this  type  has 
been  made  on  the  basis  of  a  plausible  initial  set  of  key  para- 
meters.  This  result  is  summarized  in  a  later  section. 
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6.   Summary  of  Principal  Calculation  Formulas 

The  principal  equations  developed  in  the  preceding  sections 
are  summarized  below  for  easy  reference.   They  are  listed  without 
further  explanation  in  the  approximate  order  in  which  they  would 
be  used  in  a  typical  calculation. 

Input  Parameters 


a  = 


'C_o   V 
vCb/\Po, 


x±2 


(Y-l) 


(6.1) 


S  = 


'£CB  Po 


(6.2) 


i>   =  - 


111 

Y-l 


Y+3 


(1+a) 


(6.3) 


v  =  y^3 


(6.4) 


V. 


■±-   =  (hQ-h1)   (where  s±   =  sQ,  p±   =  p4 


(5.5) 


*-  h    ,  (W  (Sc"S4) 
Ke     i1       (h  -hj 


o   1' 


(6.6) 


K 


(Pq-P4)V4 
(ho-hx) 


(6.7) 


kq    =  ivy        *, 

(h  -hn)    Kp 
o   1 


(6.8) 
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[h  -h.)    K  *, 

(no-hx)   np  f 

n  *  =   ~*  /k  *-k  *\  (6.10) 

V  I  e    p 


Solution  of  Energy  Balance  Equation 

A   =  2(y-l)2 

Al      nv  (6.11) 

A2  -fnj-  i)y2  (6-12) 

A3  =   vy2  (6.13) 

A.  =  /c. +  1)  (6.14) 

4      f    nE     / 


B4  =  -4A1-4 (3-2y)A2-5A3-4A4  (6.15) 

B3  =  4(y-l)A1-4y(y-2)A2+(10y-l)A3+4(2y-l)A4  (6.16) 

B2  =  (4y-l)A1-(4y2-2y-5)A2-y(5y-2)A3-(4y2-8y+l)A4  (6.17) 

Bl  =  YA1-(y2+l)A2~y2A3"2y(2^*1)A4  (6.18) 

BQ  =  -A2-y2A4  (6.19) 


G   =   B4x  +B3x3+B2x  +B,x+B  (6.20) 

G»  =   4B4x3+3B3x2+2B2X+B1  (6.21) 

Gn 
x  +1  =   xn  ~  G- '      (Iterate  until   x   converges)        (6.22) 
n 
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Output  Parameters 


K 


K         =  ,n  v  il±5xl   _  _P_ 

'  X   (l+2x)2    % 


K 


K 


xfy-x)    x^ 


cv 


K 
i 

k" 


PR 


Net  Electric 
Power 


Gross  Electric 
Power 


Ideal  Pressure 
Rise 

Conversion 

Efficiency 


(6.23) 


Cycle  Efficiency   (6.24) 


Relative  Output    (6.25) 


(6.26) 


(6.27) 


(6.28) 


Note:   Power  coerficients  denoted  by  symbol  K   are  expressed 
in  terms  of  the  kinetic  power  at  station  1.   Those  denoted  by 
symbol   c   are  expressed  in  terms  of  the  kinetic  power  at  station 
3.   The  conversion  between  these  two  forms  is  simply 


K   = 


(6.29) 


xy 
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7.   Some  Typical  Experimental  Data  and  Preliminary  Numerical 
Results 


Reference  (1)  reports  measurements  of  breakdown  field 

strength  as  a  function  of  pressure  at  constant  temperature  and 

5 
electrode  separation.   Their  results  show  that   C   =  9  x  10   v/m 

and  C_  =  9 . 5  x  10   m   K/C  .   Note  the  interesting  fact  that 
steam  and  air  appear  to  have  identical  breakdown  properties! 

A  system  now  under  development  by  Marks  Polarized  Corporation 
is  intended  to  operate  with  steam  at  about  100  atmospheres  pres- 
sure and  sufficient  superheat  to  give  about  5%  moisture  at  state 
1.   The  data  in  Table  7.1  suggest  that  under  these  conditions, 
approximately,  from  Eqs .  (6.1)  and  (6.2), 

1 


a   =! 


Y+3 


(Y-l) 


0.920    x   10° 

591 

JL.273] 

0.953    x   104"1 

10. 

2- 

I    x   10* 

=  0.0070  a  0 


(7.1) 


6  = 


ZCB    Po 


-B 


8.85  x  10 


-12 


k. 


n  2  r 


953  x  10 


10.34  x  10 


!J 


[591] 


=   0.0238 


(7.2) 
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Then  from  Eq.  (6.3)  we  obtain 


Y  +  l 
Y-l 


2.  3> 

0.3 


Y+3 
32 


2     3 
(1+ar  = 


1.3 


4.3 


(1+0) 


(7.3) 


In  the  absence  of  adequate  test  date,  the  best  estimate  of 
]j   that  is  possible  at  this  time  is  simply 


(7.4) 


According  to  Eq.  (6.4)  this  gives  the  estimate 


V  =  uip6  =  (1)  (1.32)  (0.0238)  =  0.032 


(7.5) 


The  value  of  the  flow  loss  factor   cf   for  the  annulus  may 
be  estimated  from  the  data  in  Ref.  (2)  on  losses  in  pipe  bends. 
At  a  Reynolds  number  of  200,000  or  above,  these  losses  are  mainly 
a  function  of   R/r,   where   R   is  bend  radius  and   r   is  pipe 
radius.   For  the  torus  in  Marks'  design   R/r  =  5.7.   Ref.  (2) 
shows  that  at  this  value  of   R/r   the  loss  factor  for  a  180 
pipe  bend  is  about  0.4.   Thus  for  a  360   turn  the  loss  factor 
would  presumably  be  twice  the  value  or  approximately 


c^  *  0  . 


(7.6) 


There  are  qualitative  grounds  for  hoping  that  the  loss  in 
a  true  torus  might  be  somewhat  smaller  than  the  above  pipe  bend 
data  would  suggest,  but  there  is  no  hard  data  available  on  this 
point.   Thus  Eq.  (7.6)  represents  a  plausible  and  conservative 
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working  hypothesis  that  can  be  advanced  at  this  time. 

Values  of   nE   can  be  estimated  from  available  test  data 
on  ejector  performance.   See,  for  example,  Ref.  (3).   These 
data  have  not  yet  been  reviewed  critically.   All  we  can  say  at 
present  on  the  basis  of  general  engineering  judgment  and  experi- 
ence is  that  for  the  range  of  ejector  design  parameters  of  in- 
terest in  the  present  context,   nE   is  unlikely  to  lie  above 
0.9  or  below  0.5.   It  is  thought  reasonable  as  a  tentative  ini- 
tial working  hypothesis  to  choose 

nE  «  0.8  (7.7) 

Also  required  for  the  calculation  are  certain  thermodynamic 
properties.   These  are  summarized  in  Table  7.1.   These  data  hap- 
pen to  be  in  English  units  Out  this  is  permissible  because  the 
final  results,  being  dimensionless ,  are  independent  of  the  system 
of  units  used.   For  some  purposes  the  steam  was  treated  approxi- 
mately as  a  perfect  gas  with  a  ratio  of  specific  heats  y   =  1.3. 

The  information  in  Table  7.1,  along  with  the  formulas  of 
Section  6,  permit  the  various  key  input  parameters  to  be  esti- 
mated as  summarized  in  Table  7.2.   These  values  are  regarded  as 
the  best  estimates  that  are  possible  at  the  present  time.   These 
input  data  suffice  to  permit  a  corresponding  performance  calcula- 
tion to  be  carried  out.   This  has  been  done  and  the  main  results 
are  summarized  in  Table  7.3  and  in  Fig.  7.1.   Notice  that  each 
row  or  Table  7.3  corresponds  to  a  different  value  of  the  ejector 
area  ratio    A, /A  . 
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Table  7.1   Thermodynamic  Properties 


State  0         14      Ambient 


p 

psia 

1500 

1100 

1100 

T 

°R 

1063.9 

1016 

1016 

h 

Btu/lbm 

1180.7 

1156.3 

557.4 

s 

Btu/lbm°R 

1.3470 

1.3470 

0.7575 

P 

lbm/ft3 

— 

— — 

45.45 

530 


Table  7.2   Estimated  Values  of  Key  Input  Parameters 


Y   =   1.3  np  =  0.9 

U   =   1  Ke*  =  12.7 

v   =   0.032  K    =  0.067 


P 


* 


HE  =   0.8  KQ   =  25.5 

cf  =   0.8  K    =  25.5 

nx  =1  nc*  =  0.50 
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Fig.  7.1   Results  of  Energy  Balance 
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The  curve  ABCDEF  shown  in  Fig.  7.1  is  the  locus  of  points 
x,y  which  satisty  the  basic  energy  balance  relation  of  the 
ejector  and  fluid  flywheel,  that  is,  Eq.  (4.23).   In  view  of 
the  tentative  nature  of  the  values  assigned  to  parameters   v, 
nE   and   c^,  the  accuracy  of  this  curve  is  somewhat  uncertain, 
but  its  general  qualitative  features  are  undoubtedly  correct. 
Any  straight  line  through  the  origin,  such  as  line  OEB,   for 
example,  represents  Eq.  (4.25),  and  the  inverse  of  its  slope 
represents  the  area  ratio  at/A3  •   Notice  that  such  a  straight 
line  may  in  general  intersect  the  curve  at  two  points  such  as 
E   and   B  ,  for  example.   It  can  be  shown  that  only  the  lower 
point  E   represents  stable  equilibrium;  point   B   represents 
unstable  equilibrium  and  should  be  disregarded.   For  this  reason 
the  entire  dashed  portion  ABC   of  the  curve  must  also  be  dis- 
carded.  Only  the  solid  portion   CDEF  has  any  practical  signifi- 
cance . 

Notice  that  the  tangent  point   C   represents  the  lowest  value 
of  A.. /A_   and  the  highest  value  of   y   that  can  be  achieved  in 

steady  state  operation.   Recall  that  the  gross  electric  power  is 

2 

proportional  to  the  factor   (l+5x) / (l+2x)  .   It  can  be  shown  that 

this  factor  reaches  its  peak  value  of  1.04  at  x  =  0.1  .   It  is 
therefore  desirable  to  choose  the  operating  point  on  the  curve 
as  close  to  x  =  0.1   as  possible.   In  this  case  the  nearest  point 
is  at   D  where  x  =  0.36.   Note  that  this  is  the  lowest  value  of 
x   that  can  be  attained  for  the  assumed  values  of   v ,  ti_   and  cf 
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The  corresponding  value  of  the  factor   (l+5x) / (l+2x)    is  0.95. 
Notice  from  the  data  in  Table  7.2  that  an  ideal  reversible 
device  operating  on  the  given  basic  thermodynamic  cycle  would 
produce  a  dimensionless  net  power  output  of  amount 

(1^*  -  Kp*)  =  12.7  -  0.07  «  12.6  (?  Q) 

■ 

and  would  produce  an  ideal  cycle  efficiency  of  amount 

nc*  =  0.50  (7.9) 

These  theoretical  performance  limits  of  the  ideal  reversible 
device  are  in  themselves  very  favorable.   Unfortunately,  the  per- 
formance of  the  actual  EHD  system  as  summarized  in  Table  7.3  is 
seen  to  be  very  unsatisfactory.   In  fact  the  net  power  outputs 
and  cycle  efficiencies  are  negative  for  all  possible  operating 
points!   This  means  that  even  at  100  atmospheres  pressure,  the 
phenomenon  of  electrical  breakdown  so  severely  limits  the  gross 
electrical  power  output  that  it  it  insufficient  even  to  meet  the 
demand  for  power  to  drive  the  pump,  let  alone  provide  any  useful 
output. 

It  might  be  supposed  that  a  radical  increase  in  the  electrical 
breakdown  strength,  if  it  could  be  achieved,  would  in  itself  trans- 
form this  pessimistic  picture  into  an  optimistic  one,  but  this  is 
not  entirely  the  case,  at  least  not  for  the  present  EHD  scheme. 
Suppose  for  example,  that  the  breakdown  parameter   v   could  be  in- 
creased by  a  factor  of  roughly  30  so  as  to  bring  the  gross  elec- 
trical power  coefficient  up  say  to 

Ke  «  1  (7.  10) 
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In  that  case  the  cycle  efficiency  would  still  be  only  about 


* 
nc  ■  £j  (Ve  "  ^  ) *  2575   (1  "  0.07)  *  0.037         (7.11) 

or  less  than  4%.   The  main  reason  for  this  unsatisfactory  result 
is  that  the  present  cycle  squanders  available  energy  by  dumping 
heat  from  the  condenser  at  high  temperature  to  the  ambient  air 
at  low  temperature  without  extracting  any  useful  power  from  it. 
The  irreversibilities  inherent  in  the  ejector  also  waste  available 
energy.   Moreover,  friction  losses  in  the  fluid  flywheel  appear 
to  be  disappointingly  high.   As  a  reflection  of  these  losses, 
note  from  Table  7.3  that  the  conversion  efficiency   n     remains 
below  4%  for  all  possible  stable  operating  points. 

It  is  also  worthwhile  at  this  point  to  review  the  various 
factors  that  affect  the  important  parameter  x>      which  so  severely 
limits  the  electrical  output.   These  factors  are  revealed  clearly 
in  Eqs .  (6.2),  (6.3)  and  (6.4).   It  is  convenient  to  combine  them 
here  in  the  following  way,  assuming  a  ~    0.   Thus 


(7. 12) 


As   y   varies  from  1  to  1.67,  the  function  of  y      in  the 
curly  brackets  above  varies  from  1.82  to  0.97.   Hence  there  is 
no  possibility  of  radically  changing  the  order  of  magnitude  of 
v   by  manipulating  y    . 
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Similar  considerations  apply  to  the  quantity   £  .   It  has 
been  suggested,  see  Re.  4,  that   £   might  be  increased  by  using 
a  polarizable  aerosol.   But  the  overall  increase  that  can  be 
achieved  in  this  way  is  expected  to  be  less  than  a  factor  of  ten 
and,  as  we  have  seen,  this  in  itself  is  not  enough  to  overcome 
fully  the  limitations  of  the  present  EHD  scheme. 

An  obvious  way  to  increase   v   is  to  increase   pQ  .   We  have 
noted,  however,  that  even  at  100  atmospheres  pressure,   v   is 
still  much  too  small.   Any  increase  beyond  this  value  is  not 
feasible  for  two  reasons.   Firstly,  structural  and  practical 
problems  become  troublesome  at  such  high  pressures.   Secondly 
the  linear  relation  between  electrical  breakdown  strength  and 
fluid  density  as  given  by  Eq.  (2.5)  breaks  down  near  pressures 
of  100  atmospheres.   See  Ref.  1. 

The  two  remaining  factors  at  our  disposal  are   C_   and  TQ  . 

2 

Note  that   v   is  proportional  to  the  quadractic  factor   (C  /T_) 

The  breakdown  constant  C   is  a  property  of  the  medium.   The 
question  whether  any  aerosol  mixture  exists  that  has  a  substan- 
tially higher  value  of   CB   has  not  been  investigated.   It  is 
recommended  that  such  an  investigation  be  carried  out. 

Finally  consider  the  factor   TQ  .   In  the  present  scheme, 
charged  particles  of  low  mobility  are  created  by  condensation. 
This  requires  that   TQ   be  high  enough  to  produce  about  5% 
liquid  at  state  1.   This  restriction  on   TQ   could  be  eliminated 
by  abandoning  the  condensation  method  and  utilizing  instead  the 
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direct  injection  of  liquid  droplets  or  solid  particles.   This 

change  also  entails  changing  the  carrier  fluid  from  steam  to 

air  or  some  other  noncondensing  gas.   Note  that  a  decrease  in 

T    from  1016   R  to  530  °R  would  increase   v   by  a  factor  of 
o 

about  3.7. 

It  seems  possible  that  some  combination  of  the  above  methods 
could  perhaps  increase   v   by  a  factor  of  say  20  or  so.   We  have 
seen  that  this  in  itself  would  not  make  the  present  fluid  fly- 
wheel configuration  successful.   Nevertheless,  an  improvement  of 
this  order  might  very  well  suffice  to  make  some  alternative  EHD 
scheme  practical. 
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8.   Conclusions  and  Recommendations. 

Broadly  speaking,  one  of  the  inherent  limitations  of  EHD 
power  generation  is  the  severe  restriction  imposed  by  the 
phenomenon  of  electrical  breakdown.   Because  of  this  effect, 
the  electrical  work  that  can  be  extracted  per  unit  mass  of  work- 
ing fluid  is  very  small.   In  general  there  are  two  distinct  ap- 
proaches for  attempting  to  overcome  this  limitation.   One  method 
is  based  on  multi-staging,  that  is,  on  passing  the  same  working 
fluid  through  a  whole  series  of  successive  electrical  conversion 
sections.   The  other  method  is  to  adopt  an  ejector  system  which 
utilizes  a  primary  stream  of  high  velocity  but  low  mass  flow 
rate  to  drive  a  secondary  stream  of  low  velocity  but  high  mass 
flow  rate.   The  conversion  of  flow  energy  to  electric  work  is 
accomplished  in  the  resulting  low  velocity  stream.   Although 
electrical  work  per  unit  mass  of  low  velocity  fluid  remains  low, 
the  electrical  work  per  unit  mass  of  high  velocity  primary 
fluid  may  be  boosted  to  a  sufficiently  high  value.   The  effective- 
ness of  this  scheme  is  limited  to  some  extent  by  the  character- 
istic ejector  mixing  losses. 

In  this  report,  a  proposed  type  of  single  stage  EHD  genera- 
tor is  analyzed  which  utilized  the  ejector  principle.   However, 
this  design  combines  the  ejector  and  the  electrical  conversion 
section  into  a  single  unit;  in  particular,  charged  particles  are 
introduced  into  the  primary  stream  only  and  are  subsequently 
collected  at  the  ejector  exit  section. 

The  analysis  shows  that,  unfortunately,  this  type  of  design 
gives  extremely  poor  performance.   In  fact  the  electrical  power 
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generated  is  insufficient  to  drive  the  feedwater  pump,  let  alone 
produce  any  net  useful  output-   This  poor  result  can  be  attribu- 
ted to  the  fact  that  the  charged  particles  are  introduced  into 
the  primary  stream,  thus  generating  a  very  intense  local 
electrical  field  at  this  point. 

An  obvious  remedy  for  this  difficulty  is  to  place  the  elec- 
trical conversion  section  after  the  ejector  where  the  electrical 
field  can  be  much  lower.   The  Appendix  evaluates  this  revised 
design.   It  is  shown  in  the  Appendix  that  under  the  most  opti- 
mistic assumptions,  the  cycle  efficiency  is  less  than  4%  which 
may  be  insufficient  to  compete  with  alternative  power  systems. 

It  may  be  possible  to  multi-stage  the  fluid  flywheel 
generators  and  thereby  show  an  order  of  magnitude  improvement 
in  performance.   It  is  recommended  that  systematic  parametric 
studies  be  carried  out  to  determine  the  performance  possi- 
bilities of  such  a  multi-stage  system  and  to  estimate  optimum 
design  values  of  the  key  parameters.   It  would  also  be  ad- 
visable in  the  interest  of  orderly  development  to  restrict  this 
proposed  next  step  to  the  same  basic  medium,  namely  steam,  and 
to  the  same  basic  Rankine  cycle  as  used  in  the  previous  studies. 
Of  course  there  is  no  assurance  that  a  multi-stage  design  will 
necessarily  prove  successful  because,  even  if  the  theoretical 
cycle  performance  is  favorable,  practical  considerations  of 
bulk  and  complexity  might  nullify  the  theoretical  advantages. 
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APPENDIX:   Performance  Characteristics  of  Revised  EHD  Generator 

The  analysis  in  the  main  text  clearly  indicates  that  intro- 
ducing charges  in  the  primary  stream  only,  limits  the  thermody- 
namic state  1  to  a  relatively  high  fixed  pressure  ratio 

(p,/p   =  0.731  for  M1  =  0.71).   An  obvious  improvement  is  to 
uncouple  the  ejector  from  the  electrical  conversion  section. 
This  is  done  in  the  following  analysis  in  which  the  electrical 
conversion  section  is  located  after  the  ejector. 

The  work  in  this  Appendix  is  based  on  a  fixed  upper  pressure 
p   =  1500  psia  and  a  fixed  steam  quality  of  0.950  at  station  1 

(See  Fig.  5.1).   However,  the  back  pressure  p,  is  allowed  to  take 
on  a  large  number  of  possible  values  ranging  from  1400  psia 
down  to  14.7  psia.   At  any  given  back  pressure,  the  area  ratio 
A3/A,  is  allowed  to  vary  as  required  to  obtain  best  performance. 
It  is  found  that  best  performance  corresponds  to  the  minimum 
possible  value  of  the  mass  flow  ratio  (m , /itu)  =  x.   This  is 
governed  by  the  energy  balance  equation  for  the  fluid  flywheel. 

Three  cases  are  considered,  as  listed  in  more  detail  on 
page  2A.   Case  A  represents  the  present  best  estimate  of  the 
three  key  parameters  n/R,  cf  and  C„.   Case  B  represents  a  much 
more  optimistic  estimate  of  parameters  nR  and  cf,  perhaps  values 
which  can  never  be  achieved  practically.   Case  C  retains  these 
optimistic  values  of  nR  and  Cf  an<^  assumes  in  addition  that  the 
electrical  breakdown  constant  CR  can  be  increased  by  a  factor 
of  /TO.   Thus,  case  C  represents  optimism  plus  a  fundamental 
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REVISED  EHD  ANALYSIS 


BASIS 

*  Design  change:   electrical  power  section  now  placed  after 
ejector. 

*  Rankine  cycle  with  - 

p   =  1500  psia 
ro        r 

steam  quality  at  state  1  =  0.95 

OPTIMIZING  PARAMETERS 

*  Pressure  p,  in  fluid  flywheel 

*  Ejector  area  ratio  A, /A_ 
CASES  STUDIED 


Ejector  Effectiveness 
Friction  Factor 


Electrical  Breakdown 
7   Constant 
(nr  °K/Coulomb) 


nE 

0.8 

0.9              0.9 

Cf 

0.8 

0.1              0.1 

CB 

9.49xJ 

.039.49xl03  /l0x9.. 
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scientific  "break-through" . 

The  single  most  important  result  of  the  analysis  is  the 

curve  of  cycle  efficiency   n    versus  back  pressure   p,  ,  Fig.  1. 

c  j. 

This  result  confirms  that  the  proposed  design  revision  has  indeed 
improved  the  performance  over  that  of  the  original  version.   Re- 
call that  the  original  design  yielded  a  negative  cycle  efficiency'. 
The  revised  design,  case  A,  produces  positive  efficiencies  for 
values  of   p,   above  about  900  psia.   The  peak  efficiency  occurs 
approximately  at   p,  =  1400  psia.   Unfortunately,  this  peak  value 
is  still  much  too  low,  below  0.003.   As  expected,  case  B  shows  a 
marked  improvement  over  case  A,  but  nevertheless   n    still  re- 
mains below  0.01.   Even  the  extremely  optimistic  case  C  gives  a 
peak  cycle  efficiency  of  less  than  0.04. 

The  efficiency  curve  for  case  C  has  an  interesting  bimodal 
form  with  the  solution  above   p,  =  1000  psia  quite  different  from 
that  below  this  limit.   The  analysis  shows  that  the  ejector  cannot 
function  for  case  C  for  values  of   p,   above  about  1000  psia,  so 
the  solution  over  this  region  pertains  to  a  simple  EHD  conversion 
in  a  single  stream;  the  ejector  principle  cannot  be  utilized  in 
this  region.   It  is  ironic  that  the  highest  efficiency  shown, 
nearly  0.04,  occurs  at  about   p,  -  1400  psia  where  the  fluid  fly- 
wheel concept  does  not  even  apply!   A  second  maximum  cycle  effi- 
ciency of  about  0.026  occurs  at  p,  =  270  psia.   This  point  repre- 
sents normal  operation  which  utilizes  the  fluid  flywheel  principle 
to  good  advantage . 
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FIG  I A     EFFICIENCY  VERSUS  PRESSURE  AT  STATION  I. 
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The  results  show  that  the  final  cycle  efficiency  that  can 
be  acheived,  and  the  optimum  back  pressure  p,  and  ejector  area 
ratio  (A3/A, )  that  will  produce  this  optimum  depend  critically 
upon  the  three  key  parameters  nE,  cf  and  C„.   If  these  para- 
meters can  be  estimated  accurately  then  design  optimization  and 
resulting  performance  can  be  accurately  predicted  and  controlled, 
otherwise  not. 

The  results  also  show,  however,  that  even  under  the  most 
optimistic  assumptions  regarding  nE,  cf  and  C  ,  the  cycle 
efficiency  is  disappointingly  low  -  less  than  0.04!   It  is 
therefore  clear  that,  viewed  as  a  single  stage  device,  the  pro- 
posed EHD  generator,  even  with  the  noted  design  improvements, 
has  performance  limitations  which  are  just  too  severe  to  permit 
it  to  compete  successfully  with  alternative  power  systems. 

On  the  other  hand,  the  possibility  still  remains  that  a 
multi-stage  version  of  this  generator  might  show  a  further  order 
of  magnitude  improvement  in  performance.   It  is  not  clear  at 
this  time  whether  the  ejector  and  fluid  flywheel  concepts  can 
be  incorporated  to  advantage  in  such  a  multi-stage  design. 

The  remainder  of  this  Appendix  summarizes  the  analysis  of 
the  revised  scheme.   The  notation  and  approach  is  generally 
similar  to  the  main  text.   The  derivations,  formulas,  tabula- 
tions and  diagrams  are  complete  but  no   attempt  has  been  made 
to  integrate  them  into  a  smooth  flowing  text. 
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ELECTRIC  POWER  COEFFICIENTS 


P  -  V. 

p-  =  c   =  vy   where   y=  (^-) 
K3     e  3 

P 

p—  =  K   =  —  where   x  =  ( — ) 
w     ex 

Ki  m3 

Pl 

v  =  0.5  8  37^ x  .  >    (CASES  A  &  B) 

h   -  h. 
o     1 

p,  =  density,  lbm/ft 
(h   -  h, )  =  enthalpy  drop,  Btu/lbm 

Note:   Above  constant  0.5837  is  based  on  a  value  of  the  break- 
down constant  of: 

CB  =  9.49  x  103    m2°K/coulomb 

An  assumed  increase  in  C„  by  a  factor  /To"  would  increase  v  by 
a  factor  of  10.   This  would  give: 

pl 

v  =  5.837^ ^-t^n  (CASE  C) 

(hQ  -  hx) 
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ENERGY  BALANCE 


OPERATING  LIMIT 

The  fluid  flywheel  can  function  only  if  parameters  r\    ,    cf 
and  v  satisfy  the  inequality 

c.(-^  -  1)  +  v  (~  +  1  +  c-)<l  (1) 

f  nE         nE       f 

OPERATING  RANGE 

Let  C0  =  (-i  -  1  +  V)  (2) 

2     nE 

Ll    —  (3) 

nE 

C   =  (-i  +  1  +  cf)  (4) 

o    nE       f 

R2  =  C?  -  4C0C  (5) 

1      <i  o 

A  solution  exists  for  any  valve  of  y  between  the  limits  - 

C,  -  R 

Y-  =  i-kr )  (6) 

mm      2C,> 

C,  +  R 

w  -  (4ct-»  (7) 

MASS  RATIO  x 

The  solution  for  mass  ratio  x  is  found  by  the  following 
sequence  of  calculations. 

.  2(y  -  I)2  (8) 

nE 

A~  =  (—  -  l)y2  (9) 

nE 

A.  =  ( ^  +  1  +  c.  +  uy2)  (10) 

3      nE       f 
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B2  =  Al  +  (3  "  2y)A2  +  A3  (11) 

Bl  =  yAl  "  (y2  '  3)A2  +  2yA3  (12) 

Bo  =  A2  +  y%  (13) 

S2  =  B2  -  4B0B  (14) 
1      z  o 

B   -  S 

x  =  (-|g— )  (15) 

Best  performance  corresponds  to  the  minimum  possible  value 
of  x.   This  is  best  found  by  numerical  trial  and  error. 
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THERMODYNAMIC  PARAMETERS 

Let      p     =  power,  watts 

p  ,   =  kinetic  power  at  station  1,  watts 

Then     K     =  =*—   =  power  coefficient 
FK2 


K   =  {1  + 
e 


:t4  "  V(so  "  V 


<h0  -  hx) 


p  ■     nr  -  hn ) 


}     Reversible  Electric  Power 


Reversible  Pump  Power 


th.  -  v     % 


KQ  =  ( (h  -  h 


Heat  Input 


p"    n  v       k* 

■r-  c^--^J(ho-hi' 


Specific  Output  Watt  Sec/Kg 


n  v    K*   - 

r       X  Pi   1 


Cycle  Efficiency 


PRINCIPAL  PERFORMANCE  PARAMETERS 


Values  used  in  calculations- 


Excitation  Factor 


Pump  Efficiency 


n  v 
x 


X 


=  1.00 


n  =  0.90 
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THERMODYNAMIC  COEFFICIENTS 


Pi 

psia 

V 

KI 

KQ 

K* 
P 

(^ 

1400 

0.4063 

54.94 

110.1 

0.0855 

0.009590 

1300 

0.1721 

26.53 

53.09 

0.0778 

0.2071 

1200 

0.1001 

17.65 

35.30 

0.0732 

0.03241 

1100 

0.06274 

12.75 

25.47 

0.0667 

0.04680 

1000 

0.04244 

10.04 

19.97 

0.0617 

0.06214 

900 

0.02958 

8.151 

16.28 

0.0568 

0.07940 

800 

0.02073 

6.769 

13.53 

0.0521 

0.09973 

700 

0.01455 

5.716 

11.43 

0.0472 

0.1233 

600 

0.01010 

4.876 

9.762 

0.0424 

0.1513 

500 

0.006839 

4.178 

8.375 

0.0377 

0.1855 

400 

0.004430 

3.574 

7.179 

0.0329 

0.2288 

300 

0.002663 

3.037 

6.113 

0.0281 

0.2871 

200 

0.001391 

2.537 

5.124 

0.0229 

0. 3702 

100 

0.0005091 

2.016 

4.083 

0.0169 

0.5211 

14.696 

0.00004396 

1.  374 

2.726 

0.0088 

1.0000 

Notes : 

Ah 

"  (ho  "  V 

AH  =  (h   -  hn  )      .  .    cnc 
o     1  p,  =  14.676 

Since  system  is  closed,  atmospheric  pressure  p.  =  14.696 
psia  has  no  particular  significance  but  serves  as  a   convenient, 
if  arbitrary  lower  pressure  limit  for  the  cycle  calculations. 
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PRINCIPAL    RESULTS 


CASE  A 

CASE  B 

CASE 

c 

pi 
psia 

nc 

e/  1 

KW  SEC 
Kg 

nc 

KW  SEC 

Kg 

nc 

p;/*i 

1400 

0.0028 

3.61 

0.0030 

3.83 

0.0360 

46.14 

1300 

0.0022 

2.86 

0.0041 

5.42 

0.0308 

41.05 

1200 

0.0013 

2.45 

0.0051 

7.06 

0.0260 

36.14 

1100 

0.0013 

1.88 

0.0062 

9.05 

0.0217 

31.40 

1000 

0.0007 

1.09 

0.0072 

10.85 

0.0183 

27.59 

900 

0.0001 

0.15 

0.0080 

12.61 

0.0186 

29.21 

800 

-0.0006 

-0.96 

0.0087 

14.26 

0.0194 

31.82 

700 

0.0092 

15.70 

0.0206 

35.22 

600 

0.0094 

16.81 

0.0220 

39.35 

500 

0.0093 

17.60 

0.0234 

44.17 

400 

0.0086 

-  17.17 

0.0248 

49.50 

300 

0.0074 

15.69 

0.0257 

54.79 

200 

0.0053 

12.14 

0.0253 

53.27 

100 

0.0018 

4.73 

0.0211 

54.57 

14.696 

-0.0067 

-8.15 

0.0050 

16.65 

13A 


1400 


Fig.    2A      v      versus    pressure   at   station    1 
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Fig.    4A      Typical      y      versus      x 
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Fig.    5A     Area    Ratio    versus      p. 
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Fig.  6A  Velocity  Ratio  versus   p1 
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Fig.  7A   Minimum  Mass  Ratio  versus   p. 
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Fig.  8A   Net  Electric  Power  versus   p. 
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ELECTRIC    POWER 


12    2       2  2    2 

P  I£CB   R    p3   A3V3         eCB   R    p3 

-£-  =   °e   =   — 1 3 =  T~ 

PK3  Ip3V3   A3  V 

vl 
Also  p-,    =    p,    and  V-,    =   — 

3  1  3        y 

2      .                      £CBV"l         £CB2r2';1 
.    .   ce  -   vy     where   v   =  -y-  =   2  <h     .   h   , 

V  -  O        1 

Reqrouping  factors  - 
2 

v  =  ( — )  — - — - — -    (dimensionless) 

2     (hQ  hx) 

In  metric  units  - 
2  2 

(£C?  R   )    =   k8.854    x    10"12    25*]  [9.  4  9    x    103   ^-^  x    461=-^-} 

2       2  vm  JL  cmb       kg   KJ 

3 
=  84.73(^L) 

V 

In  English  units  - 
2  2 

(_| )  =  [34.73^mT][9.480  x  10"4  ^p] 

x  [3.281^]3[0.4536^-]2 
L      mJ        IbmJ 

=  0.5837  5tufti 

lbm 

Hence 

pl 

v  =  0.5837,7- r—  *    (dimensionless) 

(hQ  -  hx) 
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3 
where         p,  =  density,  lbm/ft 

(h   -  h,)  =  enthalpy  drop,  Btu/lbm 
also 

P  c 

_£   =  K   =  -f=—  =    - 
PT,_     e      2    x 
Kl         xy 

SUMMARY 

P,  p,      lbm/ft 

v        0.5837     _  h    {     _      Btu/lbm 
o     1     o    1 

e    c    ,,2 
p—  =   e  =  Vy 

K3 

P 
g   _  if   -  v 

PK1     e  "  y 
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ENERGY  BALANCE 

By  a  simple  revision  of  previous  work  we  obtain  the  two  key 
equations. 

thl.)    =  (*)  (1) 

c^  ♦  <^¥  -  n  -  w  +  <^>M  -  i  -  ^iii 

(y  -  x)  'E        (y  -  x)  2 

-  (cf  +  vy2}  =  0  (2) 

Eq  (2)  can  be  reduced  to  the  form: 

AJ-^}  -  A~{x  +  (1  "  X)?}-  A-  =  0  (3) 

x  Y  '  x     2     (y  -  x)2    3 

where 

A   =  2(y  -  l)2 

1  ^  (4) 

A,  =  (-^  -  l)y2  (5) 

2  nE 


_1  .  ,  . 

'E 


A^  =  (-  -=■  +  1  +  cf  +  \>y  )  (6) 


Eq  (3)  can  be  further  reduced  to  the  form 


Bnx2  -  B,x  +  B   =  0  (7) 

2.  1       O 

32  =  A-  +     (3  -  2y)A2  +  A3  (3) 

B1  =  yAx  -  (y2  -  3)A2  +  2yA3  (9) 

BQ  =  A2  +  y2A3  (10) 

Eq  (7)  may  be  solved  by  the  quadratic  formula.   The  larger 
of  the  two  roots  is  extraneous  and  may  be  dropped. 

Let   S2  =  B  2  -  4B„B  (11) 

1       2  o 

Bl  -  S 

Then  x  =  (_^_)  (i2) 

2B2 
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For  the  limiting  case  x  =  1,  Eq  (2)  reduces  to  the  form: 
c0y2  -  C,y  +  C  =  0  (13) 


where 


■V  o 


C9  =  (-i  -  1  +  NO  (14) 

2     nE 

1    nE 


C   =  (-i  +  1  +  c-)                     (16) 

o    nE  f 

Let           R2  =  Ci2  -  4C-C                       (17) 

x  2  o 

The  corresponding  roots  are: 


(18) 

(19) 

2      2 
The  aoove  roots  are  real  only  if  R   is  not  negative. 

This  will  be  the  case  if  and  only  if: 

c.(—  -  1)  +  v  (—  +  1  +  c-)  <  1  (20) 

f  nE  nE         f   - 


min 

C,    -    R 
(    2C2 

max 

C       +    R 

(    X 
V    2C0 
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